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High energy neutrino emission from GRBs is discussed. In this paper, by using the simulation kit 
GEANT4, we calculate proton cooling efficiency including pion-multiplicity and proton-inelasticity 
in photomeson production. First, we estimate the maximum energy of accelerated protons in GRBs. 
Using the obtained results, neutrino flux from one burst and a diffuse neutrino background are 
evaluated quantitatively. We also take account of cooling processes of pion and muon, which are 
crucial for resulting neutrino spectra. We confirm the validity of analytic approximate treatments 
on GRB fiducial parameter sets, but also find that the effects of multiplicity and high-inelasticity 
can be important on both proton cooling and resulting spectra in some cases. Finally, assuming that 
the GRB rate traces the star formation rate, we obtain a diffuse neutrino background spectrum from 
GRBs for specific parameter sets. We introduce the nonthermal baryon-loading factor, rather than 
assume that GRBs are main sources of UHEGRs. We find that the obtained neutrino background 
can be comparable with the prediction of Waxman & Bahcall, although our ground in estimation is 
different from theirs. In this paper, we study on various parameters since there are many parameters 
in the model. The detection of high energy neutrinos from GRBs will be one of the strong evidences 
that protons are accelerated to very high energy in GRBs. Furthermore, the observations of a 
neutrino background has a possibility not only to test the internal shock model of GRBs but also 
to give us information about parameters in the model and whether GRBs are sources of UHEGRs 
or not. 

PACS numbers: 95.85.Ry, 98.70.Rz, 25.20.-x, 14.60.Lm, 96.50.Pw, 98.70.Sa 


I. INTRODUCTION 

Gamma-ray bursts (GRBs) are the most powerful ex¬ 
plosions in the universe. The observed isotropic energy 
can be estimated to be larger than 10®^ ergs [HQ- The 
high luminosity and the rapid time variability lead to the 
compactness problem. This problem and the hardness 
of observed photon spectra imply that 7-ray emission 
should be results of dissipation of kinetic energy of rela¬ 
tivistic expanding shells. In the standard model of GRBs, 
such a dissipation is caused by internal shocks - internal 
collisions among the shells (see reviews e.g., dii) 
Internal shock scenario requires a strong magnetic field, 
typically 10“^ — 10^ G, which can accelerate electrons to 
high energy enough to explain observed 7-ray spectra by 
synchrotron radiation. Usually, Fermi acceleration mech¬ 
anism is assumed to be working not only for electrons 
but also for protons, which can be accelerated up to a 
high energy within the fiducial GRB parameters. Phys¬ 
ical conditions allow protons accelerated to greater than 
10^° eV and energetics can explain the observed flux of 
ultra-high-energy cosmic rays (UHEGRs), assuming that 
similar energy goes into acceleration of electrons and pro¬ 
tons in the shell iHii Such protons accelerated to 
the ultra-high energy cannot avoid interacting with GRB 
photons. This photomeson production process can gen- 
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erate very high energy neutrinos and gamma rays |ir| . 
Whether GRBs are sources of UHEGRs or not, internal 
shock models predict the flux of very high energy neutri¬ 
nos at the Earth [Ullllll Many authors have investi¬ 
gated neutrino emission from GRBs. Ice Cherenkov de¬ 
tectors such as AMANDA at the South Pole have already 
been constructed and are taking data [I1II1II3. Now, 
the future Ikm^ detector, IceCube is being constructed 
[Ullllli In the Mediterranean Sea, ANTARES and 
NESTOR are under construction . If the prediction is 
correct, these detectors may detect these neutrinos cor¬ 
related with GRBs in the near future. 

In this paper, we execute the Monte Carlo simulation kit 
GEANT4 and simulate the photomeson production 
that causes proton energy loss. As a result, we can get 
meson spectrum and resulting neutrino spectrum from 
GRBs quantitatively. In our calculation, we take into 
account pion-multiplicity and proton-inelasticity, which 
are often neglected or approximated analytically in pre¬ 
vious works, although they may be important in some 
cases 1^ I23. We also take into account the synchrotron 
loss of mesons and protons. These cooling processes play 
a crucial role for the resulting neutrino spectrum. Our 
models and method of calculation are explained in Sec. 
im One of our purposes is to seek physical conditions 
allowing proton to be accelerated up to ultra-high en¬ 
ergy. Similar calculations are carried by Asano [I3, in 
which the possibility that a nucleon creates pions mul¬ 
tiple times in the dynamical time scale is included but 
multiplicity is neglected. Using obtained results, we also 
investigate efficiency of neutrino emission for various pa- 
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rameter sets. Since there are many model uncertainties 
in GRBs, such a parameter survey is meaningful. Our 
final goal is to calculate a diffuse neutrino background 
from GRBs for specific parameter sets. A unified model 
for UHECRs from GRBs is very attractive iHlIi , al¬ 
though it requires GRBs being strongly baryon-loaded 
. Even if observed UHECRs are not produced mainly 
by GRBs, high energy neutrino emission from GRBs can 
be expected. Hence, we leave the nonthermal baryon- 
loading factor as a parameter. In this paper, assuming 
GRBs trace star formation rate, we calculate a neutrino 
background from GRBs and compare our results with 
the flux obtained by Waxman & Bahcall Q. The design 
characters of neutrino detectors are being determined in 
part by the best available theoretical models, so numer¬ 
ical investigation of high energy neutrino fluxes should 
be important. Finally, we will consider the implications 
of neutrino observations and discuss a possibility that 
neutrino observation gives some information on the non¬ 
thermal baryon-loading factor and the inner engine, if 
the internal shock model is true. Our numerical results 
are shown in Sec. cni Our summary and discussions are 
described in Sec. cn 


II. METHOD OF CALCULATION 
A. Physical Conditions 

Throughout this paper, we consider the raoch that in¬ 
ternal shocks occurs (see reviews e.g., 0,0,0). We focus 
on long GRBs, whose duration is typically ^ 30 s and 
they are likely to be related with the deaths of short¬ 
lived massive stars (see Sec. ITTeIi . The internal shock 
model can produce the observed highly variable tem¬ 
poral profiles In this model, the radial time 

scale is comparable with the angular time scale (see e.g., 
0). and each time scale is related to the pulse width. 
Widths of individual pulses vary in a wide range. GRB 
pulses with ~ 0.1 — 10 s durations and separations for 
bright long bursts are typical ones . The shortest 

spikes have millisecond or even sub-millisecond widths 
in some bursts. It will reflect the intermittent nature of 
the fireball central engine. The internal collision radii 
are determined by the separation of the subshells, which 
is written by d in the comoving frame. Since the ob¬ 
served pulse width is given by St « d/2Tc, the internal 
collision radius is expressed by the commonly used rela¬ 
tion, r « 2r^cSt « 6 X 10^^(r/300)^((5t/0.1 s) cm, where 
T is a bulk Lorentz factor of the shell, which is larger 
than 100. We consider collision radii in the range of 
(10^^ — 10^®) cm in our calculation. The width of the 
subshells is constrained by the observed variability time 
and the duration, and it is typically written by I « r/T, 
where r is the radius at which internal shocks begin to 
occur. In this paper, we leave I as a parameter because 
the precise width of each subshell is unknown. The dy¬ 
namical time of each collision is t^yn ~ l/c. For sim¬ 


plicity, we assume each collision radiates similar energy, 
where N is the number of collisions and 
it is almost the number of subshells [ 23 . Although we 
do not know about N precisely, we set N ^ (10 — 100) 
in our calculation since the number of pulses per bright 
long burst is typically the order of dozens . 

The geometrically corrected GRB radiation energy is typ¬ 
ically Ery^tot = 1-24 X ergs 00, and we fix E^^tot 

throughout this paper. The isotropic energy is a few or¬ 
ders of magnitude larger than this energy, and we take 
= fb 'Ey, tot - ( 10^2 _ 1054) ergs. Here, ft = 0|/2 
is a beaming factor of GRB jet. In the prompt phase 
l/T < 9j has to be satisfied. Once we determine E^ and 
A, fb is also determined. 

The photon energy density in the comoving frame is given 

by, 


Tpiso 

u - ^ 

dTrTr^Z 


( 1 ) 


In the standard model, the prompt spectrum is explained 
by the synchrotron radiation from electrons. For simplic¬ 
ity, we assume that the GRB photon spectrum obeys a 
power-law spectrum. That is, 

dn f (for e™™ < e < s^) , . 

de (for < e < 


The observed break energy is ^ 250 keV, which cor¬ 
responds to the break energy in the comoving frame, 
~ a few keV. So we set to 1 keV for E'!y^° = 
2 X 10®^ ergs, 3 keV for Ely‘^° = 2 x 10®^ ergs, 6 keV for 
E-^so — 2 X 10®^ ergs. In addition, we set spectral indices 
to a = I, /3 = 2.2 as fiducial values and a = 0.5, /3 = 1.5 
as a flatter case. We take the minimum energy as 1 eV 
and the maximum energy as 10 MeV. Below 1 eV, the 
synchrotron self-absorption will be crucial, while above 
10MeV, the pair creation will be important [^ . 

To explain the observed emission, the internal shock 
model requires the strong magnetic energy density, which 
is expressed by Ub = ^bUj. We take cb = 0.1,1,10 as a 
parameter. 


B. Acceleration and Cooling of Proton 

The radiating particles are considered to gain their en¬ 
ergies by stochastic acceleration. The most promising ac¬ 
celeration process is Fermi acceleration mechanism. We 
assume that electrons and protons will be accelerated in 
GRBs by hrst order Fermi-acceleration by diffusive scat¬ 
tering of particles across storng shock waves. The accel¬ 
eration time scale is given by Gcc ~ /SB)^, 

where is Larmor radius of the particle. Pa = va/c, va 
is Alfven velocity, and SB is the strength of the turbu¬ 
lence in the magnetic field . Assuming that the 

diffusion coefficient is proportional to the Bohm diffusion 
coefficient, we can write the acceleration time of proton 
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as follows, 


, _ rL Sp 

= "T = ’'7S-C 


( 3 ) 


In usual cases 77 ^ 10 would be more realistic, and 77 ~ 1 
will give a reasonable lower limit for any kind of Fermi 
acceleration time scale. As for the ultrarelativistic shock 
acceleration, the acceleration time scale can be written 
tacc ^ £plV^TreieBc, where T^ei is the relative Lorentz 
factor of the upstream relative to the downstream [^ . 
For the case of mildly relativistic shocks, 77 ~ 1 is prob¬ 
ably possible [ 3 . Hence, we set 77 = 1 optimistically in 
our calculation. 

Corresponding to tacc < tdyn, proton energy is restricted 
by Ep $ eBl/rj, which satisfies the requirement that pro¬ 
ton’s Larmor radius has to be smaller than the size of 
the acceleration sites. Even when GRBs are jet-like, the 
transverse size will be larger than the radial size r/F 
because the jet opening angle 6j is larger than 1/F, so 
tacc < tdyn holds. Proton’s maximal energy is also con¬ 
strained by various cooling processes, for example, proton 
synchrotron cooling, inverse Compton (IC) cooling, and 
photohadronic cooling and so on. First, the synchrotron 
loss time scale for relativistic protons is. 


3m 1 1 
Aarml EpUs 


( 4 ) 


Second, we have to consider IC cooling, but we ignore 
IC process in this paper for simplicity. This contribution 
is considered to be small in the region where proton en¬ 
ergy is above around 10^"^eV 0 Although IC will be 
more important than synchrotron loss when proton en¬ 
ergy is below that energy, we can ignore this fact because 
photohadronic loss will dominate IC for our parameters. 
Furthermore, in our cases IC does not matter for the pur¬ 
pose to determine the proton maximum energy. Third, 
the photohadronic cooling process is important, and pi- 
ons and resulting neutrinos appear through pion decays. 
When the internal shock occurs at small radii, photo¬ 
hadronic cooling may dominate synchrotron cooling. We 
will discuss this process in the next subsection. Finally, 
we take into account the adiabatic cooling process, which 
has a time scale tad independent of the proton energy, and 
direct ejection of protons from the emission region. The 
latter may be dependent on the proton energy if diffu¬ 
sive losses are relevant. For simplicity, we neglect diffu¬ 
sive losses and assume that protons are confined over the 
time scale set by adiabatic expansion. From above, the 
total proton loss time scale is expressed as follows. 


= t-^ 

ip — ipj 


_|_ f _|_ J. i- ^ J. i- 

' ^syn ' ^ad — ^acc 


( 5 ) 


The maximum energy of proton is determined by this 
inequality, but we do not know the minimum energy. We 
set the minimum energy of proton to 10 GeV. Proton 
energy density can be expressed by Up = CaccU-y, where 
Eacc is a nonthermal baryon-loading factor. However, we 


have few knowledge about how many protons can be ac¬ 
celerated for now. So, we adopt Cacc = 10 as a fiducial 
parameter, assuming that a significant fraction of rela¬ 
tivistic protons can be accelerated. It is also assumed 
that proton distribution has dup/dEp oc according to 
first order Fermi acceleration mechanism. 


C. Photomeson Production 


GRB photons produced by synchrotron radiation and 
relativistic protons accelerated in the internal shocks can 
make mesons such as tt^ and K^, and secondary elec¬ 
trons, positrons, gamma rays, and neutrinos through the 
photomeson production process. Some authors used A- 
approximation for simplicity, but this approximation un¬ 
derestimates proton energy loss in the high energy re¬ 
gion. This is because the theoretical model predicts the 
increase of cross section in the high energy regions. We 
use the Monte Carlo simulation toolkit GEANT4 [^ . 
which includes cross sections up to 40 TeV |^. Above 
40 TeV, we extrapolate cross sections up to the maximum 
energy, which is around lOPeV. We checked the accu¬ 
racy of GEANT4 photomeson cross section within a few 
percent, comparing with PDG data [^. Moreover, we 
take into account inelasticity. Around the A-resonance, 
inelasticity is about 0.2. But above the resonance, in¬ 
elasticity increases with energy, and the value is about 
0.5 — 0.7, that is consistent with Miicke et al. |^. How¬ 
ever, it seems that the current version of GEANT4 tends 
to overestimate the amount of produced neutral pions, 
leading to underestimation of the amount of charged pi¬ 
ons by a factor of ~ 3. So we approximate 7 r+ : 7 r° = 1 : 1 
for single-pion production, and 7 r“'' 7 r“ : = 7:4 for 

double-pion production . 

Photomeson cooling time scale of relativistic protons is 
given by the following equation for isotropic photon dis¬ 
tribution M- 



where e is the photon energy in the rest frame of proton, 
7 p is the proton’s Lorentz factor, Kp is the inelasticity 
of proton, and Eth is the threshold photon energy for 
photomeson production in the rest frame of the incident 
proton, which is Eth ~ 145 MeV. Using this equation, we 
can obtain photomeson production efficiency. Protons 
with 1 PeV energy will effectively interact with soft X- 
ray photons with energy ^ O.lhkeV. 

Pion spectrum can be obtained by executing GEANT4 
through the following equation, 

max max 

dn-jr [ ^ j f d f dap-y{E,U, Ep)^ 

dEp jgTnin ^ dE J AtT dE^^ ^ ^ 

( 7 ) 

where dnp/dEp and dn/dE is proton and photon distri¬ 
bution in the comoving frame, ^ is the pion-multiplicity. 
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and tp is the proton loss time scale, which is defined by 
the equation ©. 


D. Neutrinos from Pion and Muon Decay 


By using GEANT4, we can get pion spectra. Hence, 
neutrino spectra follow from the spectra of pions and 
muons. Neutrinos are produced by the decay of —> 

+ Ve{ve) + Vp, + Vp. The life times of pi¬ 
ons and muons are tir = and tp = jpTp respectively. 
Here, = 2.6033 x 10“®s and Tp = 2.1970 x 10“®s are 
the mean life times of each particle. When pions decay 
with the spectrum, dn-j^ldeT^ by ^ + i/p{Dp), the 

spectrum of neutrinos are given by |^. 

dn„ _ m^c 1 dn-^ 

ds^, 2e* min PtT dSTT 

Here, e* = {rri^ — m^)c^/2m^ and e™ = (e*/ei/ -I- 
e;^/e*)/2. Similarly, we can get muon spectrum from 
pion spectrum. Muon decay is the three-body-decay pro¬ 
cess, which is slightly more complicated than the case of 
two-body-decay. Given the spectrum of muon, it can be 
calculated by the following equation 


dui, 

de„ 



1 dup 
cpp dsp 


W:4(/o«)Tcos0:/i(e:)) 


(9) 


where e*p^ = - ( 7 ^ - \ £*p 2 = min( 7 ^ej, -P 

( 7 ^ — 1 ) ' £y, (m^ — m1)c^/2mp), and £* is the energy 
of muon-neutrino in the muon-rest frame. fo{x) and 
fi{x) can be calculated by the quantum field theory, 
which are, fo{x) = 2x^(3 — 2x), fi{x) = 2 x ^(1 — 2 a:), 
and X = 2£l/mpC^. However, because of cooling pro¬ 
cesses of TT^ and these equations have to be applied 
at each time step. The pion or muon decay probability 
is (1 — exp(—At/t,n._^)), and meanwhile cooled by syn¬ 
chrotron cooling and adiabatic cooling. The synchrotron 
radiating time scale is given by replacing proton mass 
with pion or muon mass in the equation 0 . The adia¬ 
batic cooling time scale is still comparable to dynamical 
time scale. We neglect IG process of pions and muons in 
our calculation, because Klein-Nishina suppression will 
work in our cases 0 We also neglect neutrinos due to 
neutron decay n p + e~ + Pe, whose time scale is much 
larger than the dynamical time scale tdyn- 
Roughly speaking, in the A-resonance picture, neutrino 
spectrum follows proton spectrum above the break en¬ 
ergy which is given by, £p£^ ~ 0.3 GeV, and becomes 


harder below the break by ^ due to fp-y = tdyn/tp-y oc 

M- In addition, if > tad, neutrino spectrum 

will be suppressed by tad/tn.p due to adiabatic cooling. 
If G./j > tsyn, neutrino spectrum will be suppressed by 
tsyn/tTr,p due to synchrotron cooling. These statements 
will be confirmed numerically by above methods. 


E. GRB Diffuse Neutrino Background 

UHECRs may come from GRBs uu This hypothesis 
requires that the conversion of the initial energy of a fire¬ 
ball into UHECRs must be very efficient. In other words, 
this needs a very large value of Cacc- On the other hand, 
Waxman & Bahcall have derived a model independent 
upper bound for high energy neutrino background from 
sources optically thin to photomeson production. This 
bound is conservative and robust and it is called as 
the Waxman and Bahcall limit (WB limit). They have 
also estimated the contribution of GRBs to the high en¬ 
ergy neutrino background, assuming that UHECRs come 
from GRBs, and showed this contribution is consistently 
below the WB limit. 

Whether UHECRs come from GRBs or not, numeri¬ 
cal calculations on neutrino spectrum obtained by our 
method shown above can provide a diffuse neutrino back¬ 
ground quantitatively, assuming the GRB rate follows 
some distribution. The demonstration that long-duration 
GRBs are associated with core-collapse supernovae ^3 
implies that GRB traces the deaths of short-lived massive 
stars. Furthermore, GRBs can be detected to very high 
redshifts, unhindered by intervening dust and the current 
record is recently observed GRB 050904 This holds 
the promise of being useful tracers of star formation rate 
(SFR) in the universe 0- However, there re¬ 

main many problems such as observational bias and the 
effect of dust-enshrouded infrared starbursts, although 
these are crucial for to which extend GRBs follow the 
SFR and to which extend they can be used to determine 
the SFR at high redshifts. 

In this paper, we estimate neutrino background for sev¬ 
eral models, following Nagataki et al. |4g. First, 
we adopt the hypothesis that GRB rate traces the SFR, 
Rgrb cx Rsf- This ansatz that GRBs are likely to trace 
the observed SFR in a globally averaged sense, is not 
precluded for now although there are some uncertainties 
as described above. Second, we use parameterization of 
Porciani & Madau especially employing their mod¬ 
els SFl, SF2 and SF3, which are expressed by following 
equation in a flat Friedmann-Robertson-Walker universe. 


{z) = 0.32f^hro 


exp(3.42:) 


exp(3.8z) -I- 45 
xF{z,flrn,^A)MQ yr“^Mpc“^ (10a) 


= 0.16/*h7o 


exp(3.4z) 


exp(3.4z) -I- 22 
xF{z,Vt,ri, ^a)Mq yr“iMpc“^ (10b) 


i’* (z) = 0 . 21 /*h 7 o 


exp(3.05z — 0.4) 


exp(2.93z) -I-15 
xF{z,Vtm,^A)MQ yr“^Mpc“^ (10c) 
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where F{z, UfC) = + zf + Oa) /(I + zf^ , 

h’jQ = i?o/70 kms“^ Mpc“^, and we adopt the standard 
ACDM cosmology = 0.3, = 0.7). The correction 

factor /* is introduced for uncertainties of SFR. We set 
/* = 1 , which is consistent with mildly dust-corrected 
UV data at low redshift; on the other hand, it may un¬ 
derestimate the results of other wave band observations 
and we only have to correct /* in such cases 0. Many 
works have modeled the expected evolution of the cosmic 
SFR with redshift, but there are some uncertainties (in 
particular at high redshift z > 6 ). In addition, even ob¬ 
servational estimates at modest high redshift have been 
plagued by uncertainties arising from a result of correc¬ 
tion for dust extinction. Due to a few constraints at high 
redshift, we cannot avoid extrapolating SFRs. For these 
reasons, we employ three models. 

The inner engine of GRBs is still unknown, although 
there are several plausible candidates. One of the plausi¬ 
ble models of GRB progenitors is a collapsar model [i^ . 
because GRBs likely have a link with the explosion of 
a massive (M > (35 — 40) r otating star whose core 
collapses to form a black hole [iol 1^ 1^ . Assuming that 
massive stars with masses larger than ~ 35 Mq explode 
as GRBs, GRB rate can be estimated for a selected SFR 
by multiplying the coefficient. 


fcl X 


dm (j){m) 
Jp dm m(j){m) 


1.5 X 


( 11 ) 


where (/>(to) is the initial mass function (IMF) and mass 
is the stellar mass in solar units. Here we adopt the 
Salpeter’s IMF {4>{m) oc m“^-^®), assuming that the 
IMF does not change with time, which may be a good 
approximation if there are no significant correlations 
between the IMF and the environment in which stars 
are born. For now, extant evidences seem to argue 
against such correlations at z < 2 , although this va¬ 
lidity at high redshift is uncertain |^. For compar¬ 
ison, we can obtain i?sN(z) = 0 . 0122 Mg^' 0 *^(-z) as¬ 
suming that all stars with M > 8Mq explode as core¬ 
collapse supernovae. This result combined with /* = 1 
agrees with the observed value of local supernova rate, 
i?SN(0) = (1.2 ± 0.4) X 10-'‘h^oyr”^Mpc”^ IH. In 
the equation m, we introduce an unknown parameter, 
fcl, which expresses the fraction of the collapsars whose 
mass range is in (35 — 125)M0 are accompanied with 
GRBs. We normalize fd using the value of GRB rate, 
-Rgrb(O) = I 7 / 17 Q yr“^Gpc“^, obtained by recent analy¬ 
sis using BATSE peak flux distribution . From 

equations llOall . llObll . and llUcll . combined with m, 
we obtain following expressions in units of yr ^Gpc , 


-Rgrbi(^) = 17 


fcl 


46exp(3.4z) 

1.6 X 10-3 J exp(3.8z) -h 45 


( 12 a) 


-RgRB 2 ( 2 :) = 17 


fcl 


23exp(3.4z) 

1.6 X 10-3 J exp(3.4z) -b 22 


( 12 b) 


.RgRB3(^) = 21 


fcl \ 24 exp( 3 . 05 z- 0 . 4 ) 

1.6 X 10 - 3 / exp( 2 . 93 z) -b 15 


{ 12 c) 


We also consider the Rowan-Robinson SFR that can 
be fitted with the expression in units of yr-^Gpc”^, 


fcl 


Rrr{z) 41 ( g ^ 1 I ^qo. 75 (for 2 


10O.75Z (for 2 < 1 ) 


(13) 


The parameter, fd is unknown for now, but ~ 2 x 10-3 
will give the reasonable upper limit for the probability for 
one collapsar to generate a GRB, which corresponds to 
the value if all the GRBs come from collapsars. On the 
other hand, direct estimates from the sample of GRBs 
with determined redshifts are contaminated by observa¬ 
tional biases and are insufficient to determine the pre¬ 
cise rate and luminosity function. In addition, the ob¬ 
served peak luminosity also depends on intrinsic spec¬ 
trum. More precise data on GRBs will give us informa¬ 
tion on the accurate GRB rate. For these reasons, we 
leave the value of fd as a parameter here. 

In order to get the differential number flux of background 
neutrinos, first we compute the present number density 
of the background neutrinos per unit energy from GRBs. 
The contribution of neutrinos emitted in the interval of 
the redshift z z + dz is given as, 

dnf (E.) = Rgrb(^)(1 + zf^dz^^^^^^dElil + z)-^ 

dz 

(14) 

where = {1 + z)Ei, is the energy of neutrinos at red¬ 
shift z, which is now observed as E,^ and dNy{Ejf) /dE^^ 
is the number spectrum of neutrinos emitted by one 
GRB explosion. The differential number flux of GRB 
background neutrinos, dE,^/dE,j, using the relation 
dF,y/dEi, = cdn'ff’/dEy, 

^dN,.{{lFz)E,,) 

dEl 

(15) 

(1 -|- nmz){l F z) — Ha(2z -|- z^) 


dF,, 


dEndVt 


47ri7o 


dz i?GRB ( 
1 


where we assume z^in = 0, and Zmax = 7 or Zmax = 20. 
This is because the high redshift GRB event GRB 050904 
is recently reported by Swift observation and the GRB 
distribution likely extends beyond z = 6 • l^^ addition, 

some GRBs are expected to exist at much higher redshifts 
than z = 7. First, preliminary polarization data on the 
cosmic microwave background collected by WMAP indi¬ 
cate a high electron scattering optical depth, hinting that 
the first stellar objects in the universe should have formed 
as early as z ~ 20 |^. Second, theoretical simulations 
of the formation of the first stars similarly conclude that 
these should have formed at redshifts z ~ (15 — 40) [5^ . 
Because there are convincing evidences that at least long 
GRBs are associated with the deaths of massive stars, 
it is conceivable that high-z GRBs (z U 15 — 20 or even 
higher) exist. 
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FIG. 1: Proton cooling efficiencies, fp-y, for r = 2 x 10^^ cm 
and E'y° = 2 x 10®'^ ergs, by GEANT4 (solid line). For com¬ 
parison, we also show the case with cross section having a 
cutoff at 2 GeV (dashed line) and rough analytic approxima¬ 
tion (dotted line). 



FIG. 3: Various cooling time scales and acceleration time 
scales for case A, r = 2 x 10^® cm, and es = 0.1. The hatched 
lines show an uncertainty of the acceleration time, which cor¬ 
responds to 77 = 1 — 10 . 

III. RESULTS 

A. Photomeson Production Efficiency 

We calculate the proton cooling efficiency through pho¬ 
tomeson production by the method explained above. The 
obtained results on fp^ = tdyn/tp-y by using GEANT4 
are shown in Fig. ^ For comparison, we also show 
the case where the cross section has a cutoff at 2 GeV 
in the proton-rest frame. In this case, we checked that 
our result agrees with Asano |13| within a few percent. 
Analytic rough approximation by A-resonance is also 
shown in Fig. ^ |l(l| . In this approximation we set 
cTp.^ « 5 X 10-28 cm2, £ _ 0.3 GeV, Ae Ri 0.3 GeV, 
and Up « 0.2 in the equation ®. At 7 p ^ lO^GeV, 



log(ep [GeV]) 


FIG. 2: Proton cooling efficiencies, fp-y-, case A for r = 2 x 
10^® cm and = 2 x 10®^ ergs (solid line), case B for r = 
5.4 X 10^^ cm and E‘‘y° = 2 x 10®^ ergs (dotted line), case G 
for r = 2 X 10*^® cm and Ell‘° = 2 x 10®® ergs (dashed line). 
Each shell width is given by 1 = r/F. 



FIG. 4: Various cooling time scales and acceleration time 
scales for case A, r = 2 x 10®® cm, and cb = 1. The hatched 
lines show an uncertainty of the acceleration time, which cor¬ 
responds to 77 = 1 — 10 . 

A-resonance is a good approximation and the break en¬ 
ergy is determined by ~ 0.3GeV2. Protons be¬ 
low the break energy mainly produce pions with photons 
whose energies are above the break energy, so it leads to 
fyr oc Ep-®. On the other hand, protons above the break 
energy mainly interact with harder photons, which leads 
to /tt oc e“-®. The photomeson production efficiency ob¬ 
tained by GEANT4 is larger than other two cases which 
have the cutoff and monotonically increasing in the high 
energy region of the order of jp ^ 10^’GeV. This is be¬ 
cause Kp 0.5 —0.7 rather than Kp « 0.2 at A-resonance 
are satisfied and multi-pion production occurs in this re¬ 
gion. 

In this paper we adopt three parameter sets for GRB 
isotropic energy, photon break energy, and photon spec- 
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FIG. 5: The same as Fig. But a result for case B, r = 
5.4 X 10^"^ cm, and es = 1. 

tral indices. That is, case A: = 2 x 10®^ ergs, 

= 1 keV, and a = 1, /3 = 2.2; case B; = 

2 X 10®^ ergs, = 3keV, and a = 1, /3 = 2.2; case C: 
El^^° = 2 X 10^3gj.gg^ ^ g]^gy^ ^ ^ 0 5^ p ^ ^5^ 

Fig. 12 shows the proton cooling efhciencies for each case. 
Case A and case B have the fiducial spectral indices of 
GRB photons, while case C has a rather flatter photon 
spectrum. In addition, case C corresponds to the case 
where a more energetic burst will be observed. 

Fig. 12- Fig. El show various cooling time scales and 
acceleration time scale. Proton’s maximum energy is de¬ 
termined by equation ©. Note that most of protons will 
be depleted in the case optically thick to photomeson 
production, even if a fraction of protons are accelerated 
up to the maximal energy indicated by the equation ©. 
For example. Fig. Eland Fig. EJcorresponds to such cases. 
Only in the optically thin case, a significant fraction of 
protons with the maximum energy can escape from the 
source. In these figures, for comparison, we also show rj 
in the equation © as a parameter having the range of 1 - 
10 (hatched lines). Fig. Elshows the case where the pho- 
tohadronic cooling is crucial for determining the maxi¬ 
mum energy. Multiplicity and inelasticity will be impor¬ 
tant because the contribution of photohadronic cooling 
is comparable to synchrotron cooling if cross section has 
the cutoff at 2 GeV, or smaller than synchrotron cool¬ 
ing in the case of analytic A-resonance approximation. 
At inner radii, photon density is large, so photohadronic 
process can be important unless es is enough large for 
synchrotron cooling to be a dominant process. When the 
proton’s maximum energy is determined by the photo¬ 
hadronic process, a significant fraction of protons with 
the very high energy cannot escape from the accelerating 
site in many cases of GRB parameters. It is only pos¬ 
sible within the limited range of radii or for the case of 
the moderately smaller radiation energy than our cases. 
In many cases of possible parameter sets for GRBs, the 
synchrotron cooling is the most dominant process to de¬ 
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FIG. 6: The same as Fig. El But a result for case B, r = 
1.6 X cm, and eb = 1. 

termine the maximum energy and such a case is demon¬ 
strated by Fig. 0 and Fig. El At outer radii, the dy¬ 
namical time scale will be more and more important due 
to decreasing of photon energy density and magnetic en¬ 
ergy density. Fig. El corresponds to such a case and the 
maximum energy of proton is restricted by tad- In our 
cases, significant acceleration of protons is possible only 
at larger radii, r ^ 10^^ cm. This result is consistent 
with Asano [l3. Smaller and larger r are favorable 
to generate UHECRs. Fig. El and Fig. El demonstrate 
such cases where sources are optically thin to photomeson 
production and the production of UHEGRs is possible. 
The effects of multiplicity and high-inelasticity appear in 
the very high energy region and enhances photohadronic 
cooling of protons. However, these effects affect cooling 
time scale only for the cases of inner collision radii and 
smaller magnetic fields, cs = O.I. In the next subsection, 
we will see these effects can become important only for 
limited cases. 


B. Neutrino Spectrum and Flux 

We can get the spectrum of pions from photomeson 
production efficiencies calculated by GEANT4. From 
pion obtained spectra we can also calculate neutrino 
spectra following from the method explained in Sec. irm 
and imi Fig. 0 shows spectra of -I- Vi, as one of 
our results for r = 2 x 10^^ cm, F = 100, and F = 300. 
The high-energy break is mainly determined by tsyn/G- 

1/2 

Hence the high-energy break changes satisfying e® oc 
(see Fig. [TJ. Note, from the equation ©, different values 
of F give the similar neutrino spectra for the same E'!y^° 
and r as long as I — r/T is hold. 

Fig. E- Fig. m show spectra from single-pion, double- 
pion, and multi-pion production origins. As seen in Fig. 
El and Fig. El for the case of a = 1 and j3 = 2.2 that 
is typical for GRB, the effects of double- and multi- 
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FIG. 7: Neutrino spectra in the comoving frame with ts 
changing for r = 2 x 10^® cm. The upper lines are for F = 100 
and Elj‘° = 2 x 10®^ ergs. The lower lines are for F = 300 and 
= 2 X 10®^ ergs. Three cases are shown, es = 0.1 (solid 
line), es = 1 (dashed line), es = 10 (dotted line). 



FIG. 9: The same as Fig. |H1 But a result for r = 5.4x lO’^'^ cm 
on case B. es is set to 1. 

pion are negligible or comparable to that of single-pion, 
whose contribution can be well described by A-resonance 
approximation. Protons of energy Sp ^ 10® GeV can 
interact only with the steep part of the photon spec¬ 
trum above the break, dn/de oc Hence the con¬ 

tribution of non-resonance, SpS » 0.16 GeV, is negligi¬ 
ble. Protons with the energy Sp ^ 10® GeV can inter¬ 
act with the flatter part of the photon spectrum below 
the break, dn/de oc £~^. In this case, the contribution 
of double- and multi-pion production can be important 
because the flatter part of photon spectrum cover sig¬ 
nificant energy range. Such contribution can be crucial 
at the very high energy range, £p > 10^ GeV. At in¬ 
ner radii the contributions of double- and multi-pion are 
negligible or comparable (see Fig. |H1). Even when it 
is comparable, such a region is above the high-energy 
break. Because the nonthermal proton’s maximum en¬ 
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FIG. 8: Muon-neutrino (r'^ -I- Up,) spectra in the comoving 
frame. Muon-neutrinos are produced by the decay of pion 
and muon which origins are single- or double- or multi-pion 
photomeson production. Each case is shown. A result for 
r = 2 X 10^® cm on case A. es is set to 1. 



FIG. 10: The same as Fig. |H1 But a result for r = 1.8 x 
10®"* cm on case G. es is set to 0.1. 

ergy is around (10® — 10®’®) GeV, there are only a few 
very high energy protons which can produce multi-pions. 
At outer radii these contributions are comparable to or 
larger than single-pion fraction by a factor of ^ (2 — 3) 
(see Fig. EJ. However, in most of the region where multi- 
pion production dominates, the resulting neutrino spec¬ 
tra are suppressed by synchrotron and adiabatic cooling 
processes because such a region belongs to the high en¬ 
ergy region above around the high-energy break. For the 
case of a = 0.5 and (3 = 1.5, which is the flatter photon 
spectrum, the multi-pion production has the significant 
effect for neutrino spectra. In this case there are sufficient 
high energy photons, which can interact with very high 
energy protons, £p ^ 10^ GeV. Fig. Uni demonstrates one 
of such cases. In this case the contribution of multi-pion 
origin dominates single-pion origin by one order of mag¬ 
nitude even around the high-energy break. Note that in 
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FIG. 11: Muon-neutrino -I- Vfj) spectra in the comoving 
frame for various collision radii on case A. The subshell width 
is Z = 6.7 X 10^° cm. es is set to 1 

the very high energy region, not only inelasticity but also 
multiplicity are also high. As a result, the average pion’s 
energy which can be estimated by the parent proton’s 
energy multiplied by inelasticity and divided by multi¬ 
plicity cannot be so large [^. 

The difference between neutrino spectra from pion decay 
and from muon decay is explained as follows. Neutri¬ 
nos from muons dominate those emitted directly from 
pions in the low energy region. When a pion decays as 
TT^ ^ -f the energy fraction a muon obtains 

is typically ~ ~ 0.76. Hence a direct neutrino 

from TT^ has a fraction of ~ 0.24. When a muon decays 
as /i^ —> -h Vei^e) + + Vfi, each of three species 

will carry similar energy. But a neutrino from has a 
smaller energy fraction than a direct neutrino from tt^, 
since muon has the longer life time than pion so that it 
is subject to synchrotron cooling. 

Fig. ^2 shows neutrino spectra that occur at various col¬ 
lision radii with the fixed shell width. In this case photon 
energy density chan ges with oc r“^. This result is con¬ 
sistent with Asano [l3. Fig. 1121 shows neutrino spectra 
which occurs at various collision radii with changing the 
shell width holding I « r/F. In this case photon energy 
density changes with oc oc St~^. Here, St is 

the typical variability time scale. Roughly speaking, the 
high-energy break is determined by tsyn/tn again. This 
implies the high-energy break is proportional to 
On the other hand the low-energy break is determined 
by the minimum energy of pions produced from protons. 
The dynamical time scale, tdyn is proportional to Z, while 
tpj and tsyn are proportional to r^l, so the proton cool¬ 
ing efficiencies, which is expressed by fp-y = tdyn/tpy and 
fsyn = tdyn/tsyn, are proportional to r~^ independently 
of I [T^. So, this minimum energy will depend on only 
r. Both the low-energy and high-energy break increase 
roughly proportionally to r. Fig. IIII and Fig. 1121 confirm 
these statements. 
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FIG. 12: Muon-neutrino {uy -\- Vy) spectra in the comoving 
frame for various collision radii on case A. The subshell width 
is given by I = r/F. es is set to 1 

To get the total neutrino spectrum of a single burst 
with a collision radius r and a subshell width I simply, 
we only have to multiply the number of collisions N. 
However, here we cumulate the neutrino spectra over var¬ 
ious internal collision radii, assuming each collsion emits 
similar energy. It is because internal shocks may occur 
within some range of distances, although the distribution 
of collision radii is not so clear. In the internal shock 
model, the sequence of collsions takes place. Collisions 
of the subshells with the smaller separation will occur 
at smaller radii. The subsequent collsions will occur at 
larger radii. For simplicity, we assume the number of 
such collisions is proportional to A/d, where A is the 
total width of the shell. We also assume the fixed sub¬ 
shell width through one burst. We may need to take into 
account the spreading of subshells since internal shock 
radii are roughly comparable to the spreading radii of 
subshells 0. Even so, it does not affect the resulting 
spectra so much because photon density becomes much 
smaller at larger radii. To improve our calculation, we 
will need a detailed calculation on the internal shock 
model which is beyond the scope of this paper. 

We adopt three parameter sets in this study. That is, 
set A: = 2 x 10*^^ ergs, £*’ = IkeV, a = 1, /3 = 2.2, 

and r = (10^^ — 8.1 x 10^^) cm, and I = 6.7 x 10^° cm; 
set B: = 2 x 10®^ ergs, £*' = 3keV, a = 1, fS = 2.2, 

r = (2.7 X 10^'^ - 7.3 x 10^^) cm, and Z = 1.8 x lO^^ cm; 
set C: = 2 X lO’^^ ergs, £^ = 6keV, a = 0.5, /3 = 1.5, 

r = (10^^ — 8.1 X 10^^) cm, and Z = 6.7 x 10^° cm. Set A 
demonstrates internal shocks begin at somewhat smaller 
radii, r ~ 10^^ cm. In this set, the typical variability 
time scale is St ^ 30 ms and ii N = 200, the typical jet 
angle is 6*^ ~ 0.1 rad. Set B demonstrates internal shocks 
begin at r ^ 10^^ cm. The typical variability time of this 
set is St ~ 0.3 s and if TV = 20, the typical jet angle is 
6j ~ 0.1 rad. Set C is a very energetic case which has 
the flatter photon spectrum. Three sets are shown in 
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FIG. 13: The observed muon-neutrino -|- Vp) spectra for 
one GRB burst at 2 = 1. Neutrino spectra from pion and 
muon are shown respectively on set A with N = 200, set B 
with N = 20, set G with N = 20. Set A and Set B are es = 1, 
but Set C is cfl = 0.1 

Fig. [island Fig. |T3| To evaluate observed flux from one 
burst, we set a source at z = 1. Obtained neutrino flux of 
set A and set B are comparable with Guetta et al. [T3| . 
However, such levels of neutrino flux are hardly detected 
by km^ detector such as IceCube. Only the most pow¬ 
erful bursts or nearby sources can give a realistic chance 
for detection of Vp [ij. In our sets, only set C has the 
prospect for detection by IceCube. To see this, here we 
estimate neutrino events in IceCube. We use the follow¬ 
ing fitting formula of the probability of detecting muon 
neutrinos 

where j3 = 1.35 for Ep < 10^'® GeV, while (3 = 0.55 for 
Ep > 10"^ ® GeV. Using a geometrical detector area of 
Adet = Ikm^, the numbers of muon events from muon- 
neutrinos a burst are given by, 

Np{> A., 3 ) = Adet f dEp P{Ep) (17) 

JXeV CLl^yaA 

where Ep = 10^ GeV Hence, the numbers of muon- 
neutrinos to be expected by IceCube for set C with 
N = 20 are Np = 1.9 particles. In the case of set A with 
N = 200 and set B with N = 20, we obtain Np = 0.05 
and Np = 0.004 respectively. Of course, if eacc is more 
larger, flux can be enhanced. But too large Cacc will be 
suspicious. We will discuss this later. 

Since we obtain neutrino spectra on various parameters 
above, we can calculate a diffuse neutrino background 
from GRBs for our specific parameter sets. The results 
are shown in Fig. E)- Fig. El Fig- Elis for z^ax = 7. 
Fig. Uni and Fig. |T7| are for Zmax = 20. Our adopted 
SFR models generate similar results. When Zmax is 20, 



FIG. 14: The observed muon-neutrino {i/p + Vp) spectra for 
one GRB burst at z = 1. Neutrino spectra with various es 
values are shown respectively on set A with N = 200, set 
B with A = 20, set G with N = 20. Set A and Set B are 
es = 1, but Set G is eb = 0.1 

SF3 model gives higher flux by a factor of ~ (2 — 3) be¬ 
low the low-energy break than other SFR models. This 
is because SF3 model predicts higher SFR at high red- 
shifts. Whether Zmax is 7 or 20 does not affect neutrino 
flux up to a factor. Neutrino signals from GRBs can be 
marginally detected or not by IceGube in both figures. 
In fact, when es is set to I and we use SF3 model, we 
can obtain Np = 14 particles a year for Zmax = 7 and 
Np = 17 particles a year for z^ax = 20. On the other 
hand, in the case of set B, we get Np = 1.2 particles 
a year for Zmax = 7 and Np = 1.5 particles a year for 

-^max — 20. 

Our result for set A is very similar to the prediction of 
Waxman & Bahcall |^, although our ground is different 
from theirs. In this parameter set, our calculation on 
neutrino background from GRBs also satisfy WB bound, 
although this case is optically thick to photomeson pro¬ 
duction at inner radii, r < 10^^ cm. Set B expresses neu¬ 
trino spectra for the case of larger collision radii than 
set A. In this parameter set, GRB sources are optically 
thin to photomeson production and can be sources of very 
high energy cosmic rays. It has to satisfy WB bound and 
this implies tacc can be constrained from UHEGR obser¬ 
vations if set B is fiducial for GRBs and other parame¬ 
ters are appropriate. If GRBs are sources of UHECRs, 
this set implies GRB neutrino fluxes with z-evolution, 
E'^dFp/dEpd^l ~ 10“® GeV cm^ s“^ str“^. However, this 
suggests very large baryon-loading factor, tpcc 100 if 
current GRB rate estimation is correct and our model is 
valid. If set A is more fiducial, similar arguments leads 
to even larger baryon-loading factor because UHEGRs 
can be accelerated only at large radii. So GRBs are not 
main sources of UHECRs in set A. If we adopt larger 
isotropic energy by one order in this set with N fixed, 
the flux level will increase within a factor, because larger 
isotropic energy implies smaller Oj (but for set B, it will 
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FIG. 15: Diffuse neutrino background from GRBs for Zmax = 
7 on several SFR models. R-R means the case using Rowan- 
Robinson SFR. SFl - SF3 correspond to the models of Por- 
ciani & Madau. The left lines are for set A, while the right is 
for set B. For comparison, we show WB bounds (three dashed 
lines). The upper WB bound is for z-evolution of QSOs. The 
lower is for no z-evolution. tacc = 10 and es ~ 1 are assumed. 


FIG. 16: Diffuse neutrino background from GRBs for Zmax = 
20 on several SFR models. R-R means the case using Rowan- 
Robinson SFR. SFl - SF3 correspond to the models of Por- 
ciani & Madau. The left lines are for set A, while the right is 
for set B. For comparison, we show WB bounds (three dashed 
lines). The upper WB bound is for z-evolution of QSOs. The 
lower is for no z-evolution. tacc = 10 and eb = 1 are assumed. 
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FIG. 17: Diffuse neutrino background from GRBs for ^max = 
20 with es changing. But the large baryon-loading factor, 
Eacc = 100, is assumed. The upper lines are for set A, while 
the lower lines are for set B. All lines use the SF3 model. 


increase by one order of magnitude, which is easily seen 
from the equations o, ®, and 0). Fig. im shows 
es dependence of neutrino background. But the near 
extreme case, Cacc = 100 is presumed. If it is possible, 
neutrino will be surely observed by the detector such as 
IceCube. For example, when es is set to 1 and we use 
SF3 model, we obtain iV^ = 170 particles a year for set 
A with ^inax = 20 and Nf^ = 15 particles a year for set B 
with Zniax = 20. Observations of neutrino background, if 
detected, may give us the evidence of protons being ac¬ 
celerated in GRBs, support to the internal shock model 
of GRBs, and information about these parameters inde¬ 
pendently of X/y rays from GRBs. 


Finally, we summarize the parameter dependence of our 
results as follows. The collision radii, r and the width of 
a subshell, I determine the photon energy density. The 
larger radii and width of a subshell make the result¬ 
ing neutrino emissivity smaller. The low-energy break 
of neutrino spectrum, is determined by the photon 
break energy, but it is also roughly proportional to r. 
The high-energy break is determined by the synchrotron 
(or adiabatic) cooling and satisfies e® oc We take 

eb = 0.1,1,10 (in the case of aftergrow, eb = 0.1 is 
preferred). The high-energy break is also roughly pro¬ 
portional to r^l. One of the most important parameters 
is the nonthermal baryon-loading factor, Eacc- We set 
Eacc to 10 except Fig. El and the larger Eacc can raise 
the flux level of neutrino, although too large Eacc is not 
plausible. We take ~ (10®^ — 10^'*) ergs, and more 

energetic bursts can produce more neutrinos when other 
parameters are fixed. 


IV. SUMMARY AND DISCUSSION 

In this paper we calculate proton’s photomeson cool¬ 
ing efficiency and resulting neutrino spectra from GRBs 
quantitatively. We are able to include pion-multiplicity 
and proton-inelasticity by executing GEANT4. These ef¬ 
fects of multi-pion production and high-inelasticity in the 
high energy region enhance the proton cooling efficiency 
in this region, so they help prevent protons from acceler¬ 
ating up to the ultra-high energy region in several cases. 
But in many cases, the synchrotron loss time scale and 
the dynamical time scale determine proton’s maximum 
energy. Furthermore, resulting neutrino spectra are not 
so sensitive to the nonthermal proton’s maximum energy 
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except the very high-energy region above the high-energy 
break. 

In our cases GRBs are optically thick to photomeson pro¬ 
duction at r < 10^^ cm, while at r > 10^^ cm GRBs are 
optically thin to it, so the production of UHECRs is pos¬ 
sible at larger radii. This result is consistent with Asano 
[T^ . Using the obtained proton cooling efficiency, we can 
calculate neutrino spectra. The effects from multi-pion 
production on resulting spectra are also calculated. We 
show that the contribution of multi-pion is almost negli¬ 
gible at inner radii but can be larger than that of single- 
pion at outer radii by a factor. In addition, the contri¬ 
bution of multi-pion production is somewhat sensitive to 
the proton’s maximum energy, which would be actually 
difficult to determine precisely due to uncertainty of rj. 
But such a contribution can be significant for the flatter 
photon spectrum, even though mesons lose their energy 
through these cooling process. Radiative cooling of pi¬ 
ous and muons plays a crucial role in resulting spectra. 
Neutrino spectra are suppressed above the high-energy 
break energy. If the magnetic field is strong, such sup¬ 
pression becomes large and vice versa, if the magnetic 
field is weak, such suppression becomes weak. The ob¬ 
servations of neutrino has the possibility that gives us 
some information about such a parameter of GRBs inde¬ 
pendently of gamma-ray observation. However, as shown 
in Dermer & Atoyan only the most powerful bursts, 
which are brighter than ^ 10®^ ergs or bursts at z < O.I, 
produce detectable neutrino bursts with a km^ detector 
such as IceGube. Set G of our parameter sets would be 
one of such detectable cases, but only a few neutrinos are 
expected even for the brightest bursts. For one neutrino 
burst, we adopt T = 300 and consider only the on-axis 
observations. If we observe at off-axis, we will observe 
lower energy neutrinos. 

A diffuse neutrino background is also calculated in this 
paper. In our specific parameter sets, neutrino back¬ 
ground observations by IceGube can expect a few or a few 
tens order of neutrinos per year, although it is important 
which parameter set is fiducial. Extrapolation of SFR 
to high redshifts may not be valid. Even so, our results 
would not be so much affected as we have seen above. In 
addition, GRBs may trace not SFR, but metallicity. In 
the collapsar model, the presence of a strong stellar wind 
(a consequence of high metallicity) would hinder the pro¬ 
duction of a GRB, therefore metal-poor hosts would be 
favored sites |49| . There remains large uncertanty at low 
metallicity at present, and as more bursts are followed up 
and their environments are better studied by Swift, this 
correlation will be testable. However, our results will not 
be changed so much even when we can take into account 
this. 

Throughout this paper, we set N to the range of ~ 
(10 — 100). If we change N, a neutrino signal from one 
burst will change with N. On the other hand, when 
we calculate a neutrino background, the results are not 
changed unless we fix This is because 9j varies ac¬ 
cording to the change of N since we fix and E^^tot- 


However, the results will be changed when we fix 6j and 
E^^tot- For example, on set B, we will get the lower flux 
level by one order of magnitude if we adopt N = 200 
fixing Ej^tot, 9j, and eacc- This is easily seen from the 
and the equations QJ, ©, and |7J). 
When we fix Ej^tot, dj, and Sacc more shells mean lower 
photon energy density. For set B, which is the case opti¬ 
cally thin to photomeson production, it leads to decrease 
dn^/d£,r roughly by two orders of magnitude. Hence, 
the total neutrino spectra will be lower by one order of 
magnitude. If we fix 9j, our results in which N is set to 
~ 10, would give the reasonable flux level in the opti¬ 
mistic case. 

In our parameter sets both one burst emission and a dif¬ 
fuse neutrino background give the neutrino flux we can 
barely observe by IceGube. To raise flux, GRBs require 
the larger nonthermal baryon-loading factor, eacc, which 
is difficult to estimate from microphysics at present. 
However, there are some clues and assuming too large 
Eacc will not be plausible. First, as seen in Sec. liirm 
UHECRs observations can give the upper limit to eacc- 
Furthermore, the large baryon-loading factor suggests a 
significant contribution of the accelerated protons in the 
observed hard radiation through secondaries produced in 
photomeson production. Such emission is expected to be 
observed in the multi-GeV energy range by electromag¬ 
netic cascades [Til 1313 ■ If the flux level of multi-GeV 
emission is comparable to the neutrino flux obtained in 
this paper, the flux level will be below that of X/y emis¬ 
sion and the EGRET limit. But assumed photon spectra 
will be modified by the radiation from secondaries. In 
addtion, such high energy emission may be detected by 
the near future GLAST observation |3- Although it 
is important both to compare X/y emission with such 
multi-GeV emission and to investigate whether GLAST 
can detect such multi-GeV emission or not, a detailed cal¬ 
culation for this purpose is needed and it is beyond the 
scope of this paper. If such a calculation is done, GLAST 
observation in the near future has the possibility to give 
more information about eacc- This will be the second 
clue about eacc- Third, eacc will be constrained by the 
GRB total explosion energy, which is still unknown. The 
large baryon-loading factor leads to the large explosion 
energy. For example, if the true total explosion energy is 
< 10®^ ergs, this suggests eacc ~ 100. 

On the other hand, GRBs associated with supernovae 
may imply the isotropic kinetic energy, > 10®^ ergs, 
which is larger than usual supernovae by one order of 
magnitude [3 ■ This leads to a collapsar model as a failed 
supernova in the sense that a core collapse event failed 
to form a neutron star and instead produced a black hole 
[ 4 ^. Here, we assume the internal shock scenario is cor¬ 
rect and a collapsar model is valid. If the true total explo¬ 
sion energy of a collapsar is assumed to be < 10®^ ergs, 
neutrino spectra will be limited by fcieacc ~ 10“^ be¬ 
cause fci is constrained by fd ~ 10“^ from current GRB 
rate estimations. So if observed neutrino spectra is higher 
than these values, the collapsar model cannot explain the 
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spectra alone. On the other hand, the supernova model 
predicts higher neutrino flux by a pulsar wind . If 

the observed neutrino flux is higher than expected in the 
collapsar model, the supernova model might be likely to 
exist. 

If IceCube can detect neutrinos and confirms they have 
the expected level of neutrino flux by our calculation, this 
will be one of the evidences that a significant fraction of 
protons can be accelerated and our employed internal 
shock model is valid. Of course, to estimate neutrino 
background more precisely and make our discussion jus¬ 
tified, we have to choose the most fiducial parameters for 
GRBs. The results depend on photon density, and the 
contribution from a fraction of bursts with large photon 
density might be large. So, we should take into account 
the respective distributions of parameters to execute the 
most refined calculaion. Unfortunately, many parame¬ 
ters have large uncertainty at present. For this reason, 
we calculate for a wide range of these parameters in this 
paper. The signature of GRBs may depend on z. For ex¬ 
ample, the total isotropic energy of a GRB and photon 
spectral indices may depend on z. More and more obser¬ 
vations in the near future and more refined theoretical 
models will allow our results to be improved. 

So far, we have not taken account of neutrino oscilla¬ 
tions. Since we have considered many decaying modes, 
the production ratio of high energy muon and electron 
neutrinos is not 2:1 exactly. However, the neutrinos will 
be almost equally distributed among flavors as a result 
of vacuum neutrino oscillations [Tn|. So there may be a 
possibility that tau neutrinos are detected through dou¬ 
ble bang events |^ . 

In summary, we obtain the neutrino spectrum from 
GRBs quantitatively by using GEANT4 simulation kit. 
We show that photomeson cooling process can constrain 
the proton’s maximum energy and the effects of multi- 
pion production and high-inelasticity can enhance the 
cooling efficiency. Furthermore, these effects affect the 
resulting neutrino spectra slightly and can be signifi¬ 
cant for the flatter photon spectrum. We quantitatively 


checked radiative cooling of pion and muon play a crucial 
role, which is controlled by e^. We also confirmed that 
UHECRs can be accelerated at r > 10^'^ cm. We have 
calculated not only neutrino spectra from one burst but 
also the GRB diffuse neutrino background using several 
SFR models. We have found our specific parameter sets 
give neutrino spectra comparable with the prediction of 
Waxman & Bahcall M without supposing the very large 
nonthermal baryon-loading factor, which is necessary for 
the assumption that GRBs are main sources of UHECRs. 
We have also discussed influences on neutrino spectra by 
changing parameters. Such a study is important since 
there are many parameters. If neutrino signals are de¬ 
tected by AMANDA, ANTARES, NESTOR, or IceCube, 
it will be one of the evidences that protons can be acceler¬ 
ated to very high energy in GRBs and the internal shock 
model of GRBs are plausible. Furthermore, such neu¬ 
trino observations in the near future may give us some 
information about the nonthermal baryon-loading factor 
and the inner engine of GRBs. 
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